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We have developed a microrheometer, based on optical tweezers, in which hydrodynamic coupling between 
the probe and fluid boundaries is dramatically reduced relative to existing microrheometers. Rotational 
Brownian motion of a birefringent microsphere within an angular optical trap is observed by measuring the 
polarisation shifts of transmitted light. Data gathered in this manner, in the strongly viscoelastic fluid 
Celluvisc, quantitatively agree with the results of conventional (bulk) rheometry. Our technique will 
significantly reduce the smallest sample volumes which may be reliably probed, further extending the study 
of rare, difficult to obtain or highly nonhomogeneous fluids. 

The field of microrheology, like the macroscopic study of flow and deformation, aims to characterise the 
manner in which fluids respond to applied stresses and/or strains. For soft materials this information is most 
commonly deduced by observing the motion of one or more microscopic probes embedded in the medium 1 . 
Microrheology has been applied in the study of a variety of physical systems, including polymer networks 2 " 7 , 
colloids 2 ' 8 " 10 and fluids with mesoscopic nonhomogeneities 4,11 . These systems, collectively known as complex 
fluids, are said to be viscoelastic, as their response to deformation may exhibit characteristics of both elasticity and 
viscous dissipation. Complex fluids arise frequently in biology - examples include synovial fluid, tear films, lung 
surfactants and cerebrospinal fluid - and changes in their rheology have been linked to the development and 
diagnosis of diseased states 12 " 14 . Microrheology is therefore a powerful emerging tool for biophysics and medical 
diagnostics 1115 . 

The frequency- dependent behaviour of a viscoelastic fluid is described by the complex shear modulus, G* (co). 
The complex shear stress maybe expressed as a (co) = G* (co)s (co), where s is the complex shear strain 16 . One can 
show that G' = Re{G*}, the storage modulus, describes the amount of recoverable mechanical energy stored 
within the fluid. This is broadly similar to the interpretation of Young's modulus. Conversely, the loss modulus, 
G" = Im{G*}, reflects the momentum-dependent contribution of viscous damping. 

The extreme nonhomogeneities and small volumes encountered in the study of many biological fluids ensure 
that the determination of G* by conventional rheometry is impractical. This necessitates the development of 
accurate microrheometers capable of reliably operating within ultra-small volume (~ pL) fluid samples. Such 
microrheometers must be relatively insensitive to the effects of hydrodynamic coupling between the probe 
particle(s) 10,17,18 and any fluid boundaries, potentially separated by only micrometres, as these interactions lead 
to a position-dependent decrease in the apparent mobility of the probe 19,20 . Overcoming these hurdles will allow 
continued miniaturisation of rheometry in lab-on-a-chip-style devices. 

The influences of this coupling on the centre- of-mass and angular mobilities (in a viscous fluid) are described 
by Faxen's corrections 20 . The modification of the motion of a sphere near to a fluid boundary is specified as four 
coefficients; one each for the drag coefficients corresponding to rotations and translations parallel or perpen- 
dicular to the boundary. To leading order, the translational corrections for a particle of radius a scale as (a/s) 1 , 
whilst rotational corrections vary with (a/s) 3 , with s being the distance separating the probe's centre from the 
boundary 20,21 . As such, the rotational corrections do not become significant until s ~ a, allowing rotational 
microrheology to operate in more tightly confined geometries, or more heterogeneous fluids, than its trans- 
lational counterpart 20,21 . As we will show below, this statement holds even relatively far from a boundary, despite 
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the fact that rotational microrheology has a less favourable scaling of 
the uncertainty AG* with Aa than linear microrheology. 

Previously, single-particle rotational optical microrheology has 
been restricted to measurements of viscosity 19,22 or |G*|, rather than 
G' and G" individually. Herein we report the measurement of the full 
complex shear modulus, by a single-particle rotational microrhe- 
ometer, across a large frequency range. We achieve excellent quant- 
itative agreement between conventional rheological data and those 
gathered by monitoring the angular position fluctuations of a 
birefringent microsphere undergoing biased Brownian rotational 
motion in linearly polarised optical tweezers. 

Results 

The experimental apparatus is shown in Fig. 1. Our microrheometer 
exploits two laser sources; one forms a three-dimensional optical trap 
in an inverted microscope (Olympus, 1.3 NA, 100X objective and 
matched condenser lenses) whilst another monitors the angular 
position of the trapped probe particle. The beams are coaxial and 
share a common focal spot at the position of the trap. The single 
optical trap is generated by highly-focused, linearly polarised infra- 
red (IR) light (IPG Photonics YLR-5-1064-LP,CW, 1064 nm with 
stable output at 2.0 W). An acousto-optic modulator (AOM, 
IntraAction DTD-274HA6, driven at 27 MHz) is aligned such that 
the first diffracted order is passed into the microscope, allowing fine 
control of the trapping power. 

Vaterite (CaC0 3 ), our chosen probe material, forms highly 
birefringent (An = 0.1) polycrystalline microspherulites 23 . Despite 
non-trivial internal structure, the spheres have a well-formed 
principal optic axis, which allows the exchange of spin angular 
momentum with electromagnetic fields 24 . When held in linearly 
polarised optical tweezers these particles undergo biased diffusion 
about the stable mechanical equilibrium position, with the optic axis 
in the transverse plane and the azimuthal angle <j> 0 determined by the 
polarisation vector of the trapping beam. 

Measurements of the angular position are achieved via polarisa- 
tion tracking of a low-power (~ 3 mW at focal plane) HeNe laser 
beam (JDS Uniphase 1137P, CW, 632.8 nm) which is focused to a 
relatively broad focal spot. This tracking technique is related to those 
employed by Wood et at. and Knoner et at. 25,26 . Transmission 
through a vaterite microsphere deforms the initial circularly 
polarised beam into an elliptically polarised state, in which the long 




Figure 1 | Experimental apparatus employed to monitor the directed 
rotational diffusion of a vaterite microsphere held in linearly polarised 
optical tweezers. Labels denote the infrared trapping laser (IR, P ~ 
20 mW at the trap, X = 1064 nm), circularly polarised HeNe tracking laser 
(HeNe, P < 3 mW, X = 632.8 nm), half and quarter waveplates (X/2 or XI A 
respectively), polarising beamsplitter (P), acousto-optic modulator 
(AOM), camera (CCD), visible light source (V) and polarisation-sensitive 
photodiode detectors (Dl and D2 for HeNe and IR beams respectively). 
The projection of the HeNe beam's polarisation ellipse onto a rectilinear 
basis allows the orientation of the vaterite probe particle to be recovered. 
The AOM allows control of the transmitted laser intensity. 



axis of the polarisation ellipse runs parallel to the optic axis of the 
probe. The resulting difference in the power carried by orthogonal 
(linearly polarised) beam components encodes <j>. 

This information is recovered by passing the transmitted light 
through a polarising beam splitter cube and recording the difference 
in power arriving at each output port. Assuming the photodiode 
voltage is linear in the optical power, the recorded voltage difference 
varies sinusoidally with twice the azimuthal angle, AV — {AY) oc sin 
20, because of the probe's mirror symmetry about the optic axis. To 
recover absolute angles we observe the continuous rotation of the 
probe particle under circularly polarised IR light: once the extreme 
values of Ay are known it is possible to map the data onto the interval 
[ — 1,1] and solve for 0 as a function of AV. 

Each individual probe is calibrated separately, requiring tens of 
seconds to one minute each, to ensure that variability of size and 
internal crystal structure 23 within the population is properly 
accounted for. Note that continuous rotations induce restructuring 
of polymer networks present in strongly viscoelastic fluids, altering 
their apparent rheology or pushing them into the nonlinear regime 1 : 
this potential source of data contamination has been avoided 
by collecting passive measurements of Brownian motion before 
calibrating. 

The experimentalist concerned with small samples of nonhomo- 
geneous fluids will find this system holds several advantages over 
existing microrheometers. 

One of the simplest and most effective methods of measuring fluid 
properties is by monitoring the rectilinear motions of an ensemble of 
suspended probe particles. The well-developed field of optical multi- 
probe microrheology makes extensive use of dynamic light scattering 
(DLS) 27 and diffusing wave spectroscopy (DWS) 28 , which generally 
extract data from the correlations of backscattered light. More spe- 
cialised techniques, such as localised-coherence-volume DLS (cap- 
able of measurements localised to a volume of ~ 100 fim cubed) 29 , 
are also in use. Additionally, the capacity for rapid multi-particle 
tracking afforded by modern CCD cameras has found applications 
in microrheology 1 ' 30 . 

Since the sample volume necessarily scales at least as rapidly as the 
number of probe particles (for any given technique) it is not possible 
for any of these techniques to approach the spatial resolution (the 
minimum distance over which a difference in G* may be detected) 
achieved by a single-probe microrheometer; they are ideally- suited to 
gathering large bodies of statistical data in near-homogeneous fluids. 
To increase the spatial resolution it is necessary to turn to few- and 
single -particle microrheometers. 

Optical microrheometers based on translational motion of the 
probe particle(s) may monitor the Brownian motion of a trapped 516 
or untrapped 31 probe, observe relaxation into an optical potential 6 ' 16 or 
continuously modulate the trap position 511 , amongst other possible 
schemes. All of these measurements rely on rectilinear motion, so are 
limited in spatial resolution by the strong dependence of the trans- 
lational Faxen's corrections on the distance from a fluid boundary. 

Our rotational microrheometer increases the achievable spatial 
resolution above those attained with the methods listed above 
through two mechanisms; firstly, the centre of the probe can be 
tightly constrained without degrading the signal-to-noise ratio, un- 
like translational measurements; and secondly, the rotational Faxen's 
corrections are small even within a radius of a fluid boundary. 

Existing rotational microrheometers may also be categorised as 
ensemble or single-particle devices. In the former category are mag- 
netic rheometers, which probe G* by examining the net magnetisa- 
tion of a fluid interspersed with ferromagnetic probe beads - they are 
well-suited to measurement of low- frequency behaviour 32,33 , but are 
incapable of resolving G* at high frequencies (wlln > 1 kHz) - and 
the rotational cousin of DLS, depolarised DLS 34 . 

In the latter category are single-particle rotational microrhe- 
ometers. Many of these are faced with a trade-off between angular 
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sensitivity and ease of analysis; it is possible to optically track aniso- 
tropic and asymmetric objects, such as trapped wax micro disks 3 , but 
calculating their diffusive dynamics is highly nontrivial; conversely, 
the symmetry of spheres makes them difficult to track but simpler to 
analyse. We have circumvented this issue by employing a spherical 
particle with an optical anisotropy, allowing sensitive detection of its 
orientation with readily- available polarisers and photodetectors. 

Finally, this Brownian microrheometer is able to acquire data far 
more rapidly than a similar instrument 35 in which a birefringent 
sphere was sinusoidally driven with a well-defined frequency, 
because broad-band thermal forcing eliminates the necessity of scan- 
ning over each frequency of interest serially. A combination of these 
techniques could allow one to rapidly identify interesting features in 
G* by thermal effects, then perform high-resolution scans of those 
features with the active excitation method. 

We adopt a generalised Langevin equation to model the stochastic 
evolution of the azimuthal angle of a spherical probe with moment of 
inertia I and radius a 3 . Without loss of generality we choose 0 O = 0. 



I(j) = TTH{t)- 



£(t-t')t(t')dt'-xt 



(1) 



The angular position evolves under the effects of a Hookean optical 
restoring torque with effective stiffness the zero -mean thermal 
torque t th and the fluid's hydrodynamic contributions, here repre- 
sented by the convolution of the angular velocity with the micro- 
scopic memory function £ (t) 3 . 

The existence of a linear optical restoring torque is confirmed by 
experiments, as clearly shown in Fig. 2. The distribution of occupied 
angular positions is in accordance with a thermal ensemble within a 
harmonic potential. Therefore the angular trap stiffness may be con- 
fidently determined by way of the equipartition theorem. Note that 
this calculation is independent of the fluid's rheological properties. 

X =W) (2) 

Here brackets (•••) indicate a measurement averaged over all t; 
experimentally, the measurement time must be sufficiently long that 
the mean- squared displacement approaches a constant value. 

The normalised angular position autocorrelation function 
(NAPAF), defined by Eqn. (3), contains information concerning 
the distribution of <j> (t + t) conditioned on a measurement of 
(j) (t). As such, it encapsulates the dynamical behaviour of the probe. 




-0.05 0 0.05 

Angular Position ($, rad) 

Figure 2 | The distribution of angular positions is consistent with a 
thermally-populated effective harmonic optical potential. An 

experimental histogram of vaterite orientations (circles), collected in 25% 
Celluvisc, is shown against a normalised Gaussian (solid line) of the form 
po exp [ — 0 2 /2(0 2 )] . The close agreement indicates that the optical restoring 
torque is linear in 0 within the region sampled by the particle. 



A(t) = 



m 2 ) ' 



(3) 



By analogy with Preece et al. 16 , we apply the generalised Stokes- 
Einstein relation of Cheng and Mason 3 to Eqn. (1) to find 



G*H 



Sua 3 



Ico 2 



XcoA(co) 



i+caA(ca) 



(4) 



It is important to note that A is representative of the unilateral 
Fourier transform 1 of A, defined as 



A(CO): 



A(T)e- i0JT dT. 



(5) 



Numerical evaluation of Eqn. (5) is performed using an approxi- 
mation closely related to that of Evans et al 36 . This method treats 
the experimental data as a piecewise linear function, interpolating 
between adjacent points. Times t beyond the experimentally probed 
region are integrated over by extrapolating A (t) with an appropriate 
function. 

In order to choose said function we note that, for a simple viscous 
fluid in the low Reynolds number (R) regime 16 , Eqn. (3) reduces to 



A(t) = 



exp 



Sna 3 rj 



(6) 



The majority of materials will respond in an essentially viscous man- 
ner if deformed over timescales much greater than any characteristic 
relaxation times of their substructures, allowing us to apply Eqn. (6) 
as an approximation to the long-time behaviour of the NAPAF mea- 
sured in complex fluids. 

Given a uniformly-sampled set of N points Aj spaced by At in the 



domain % 
find 



[0, t n ], we extrapolate using A, 

icoAt^ N-l 



--A N e- a eff {T ~ TN) to 



A{w) = 



1- 



w 2 At 



-A N e 



■E 



A N e 



AAe 



- ico j At _ 



(?) 



z'c0 + a e ff 



The extrapolating function introduces only a single fitting parameter, 
the effective decay constant oc e ff. 

Much of the theory presented above is predicated on the existence 
of low-^ flows. We confirm that we are probing low-^ 5 dynamics by 
comparing the experimentally determined NAPAF (water) with its 
corresponding theoretical form in this regime, Eqn. (6). As seen in 
Fig. 3, the agreement between these is excellent. This indicates that 
Eqn. (4) may be simplified by neglecting the inertial terms ( oc Ico 2 ) 25 . 

We benchmark our method against conventional techniques by 
examining four aqueous dilutions of a common dry eye treatment, 
Celluvisc (Allergan). This fluid, a carboxymethylcellulose (CMC) 
sodium solution (10 mg mL" 1 ), is designed to mimic the rheological 
properties of the fluid film coating the cornea. CMC solutions are 
quite rheologically stable, with unpreserved samples aging over a 
number of weeks 37 . 

Microrheological measurements were conducted in dilutions of 
100, 50, 25 and 0 percent Celluvisc by weight, with trapping powers 
of approximately 20-60 mW. Rheological properties were found to 
be essentially independent of the trapping power over this range, as 
noted in previous experiments 22 . Each measurement was carried out 
over three to five minutes. 

The bulk response of each fluid was probed with a parallel-plate 
strain controlled rheometer (ARES, Rheometric Scientific). With 
time-temperature superposition at approximately 8, 25 and 40°C 
we were able to achieve measurements of G* across the range w e 
[1, 200] rad s" 1 without entering the operating regime in which the 
inertia of the instrument significantly affects data. Although these 
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Lag time (t, s) 



Figure 3 | Experimental autocorrelation functions collected in samples of 
water are in close agreement with their expected forms. A decaying 
exponential function (solid line) is fitted to the experimental data (circles), 
with a single free parameter. We find the viscosity to be rj = (0.97 ± 
0.16)mPa*s, consistent with the known viscosity near to room 
temperature. The raw data were collected in a single run of five minutes 
and are presented with no averaging or filtering. 

frequencies are comparatively low, the region of overlap with our 
tweezers-based method is sufficient for critical analysis. 

Discussion 

In optimal conditions measurements of G* were repeatable and 
consistent. The largest sources of uncertainty in the moduli stem 
from the low-frequency motion of the apparatus, which pollutes 
the NAPAF, and an inability to precisely determine the probe radius. 
The former was able to be partially mitigated by high-pass filtering 
the raw time-series prior to analysis (passing co > 1 rad s" 1 ). It may 
be possible to further suppress low-frequency noise by employing 
two (mutually incoherent) trapping beams, allowing high (rectilin- 
ear) trap stiffnesses whilst maintaining a weak angular trap. 

The latter, which introduces a particle radius uncertainty, contri- 
butes a goodly portion of the final experimental uncertainty in G*. 
With our estimated uncertainty of Aa ~ 0. la, the relative uncertainty 
in 1 1 a 3 is ~ 30%; this is smaller than the systematic error associated 
with translational wall effects if the centre of the particle is within five 
radii of the boundary. This increases to s ~ 9a if the uncertainty is 
reduced to Aa ~ 0.05a, illustrating that rotational microrheometers 
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Figure 4 | A comparison of the loss modulus of aqueous dilutions of 
Celluvisc (100, 50, 25 and 0% by weight) measured by bulk (circles) and 
microrheometric (lines) techniques. The (theoretical) loss modulus of 
water (crosses) is of the form G" = cor]. All data show excellent agreement 
over a frequency range in excess of a decade. 




Frequency (co, rad/s) 

Figure 5 | The storage modulus of aqueous dilutions of Celluvisc (100, 50 
and 25% by weight), measured by conventional (circles) and 
microrheological (lines) techniques. Water (G' = 0) is omitted for clarity. 
Apart from small deviations (discussed in text) these data closely agree. 

are capable of outperforming their translational counterparts even 
relatively far from a wall. We also note that if the same probe particle 
is used to make comparative measurements, which are of immense 
interest to the biophysicist, the uncertainty in the radius does not 
affect the results. 

Data were prepared for presentation by linearly interpolating each 
of G' and G" across four independent measurements, averaging these 
and selecting a logarithmically- spaced set of representative points. 
We do not display G* (co > 10 4 Hz) due to the presence of high- 
frequency artefacts caused by our finite sampling rate 31 ' 38 . 

Fig. 4 and Fig. 5 compare our microrheological measurements 
with conventional rheometry. The data compare favourably over 
the frequency range in which they overlap, spanning at least one 
decade. Slight deviations are evident in G' measured in 25 percent 
solution, at frequencies exceeding 10 2 rad s" 1 . Although it is difficult 
to disentangle possible contributing factors, we believe that this dif- 
ference may be due to depletion interactions 2 or the breakdown of the 
continuum approximation 3 . 

We note that the high-frequency behaviour of G" is consistent with 
shear thinning, also observed in methylcellulose solutions 5 . 

Our novel microrheological scheme has proven to reliably mea- 
sure properties of complex fluids in small fluid samples. By analysing 
the rotational motion of a single probe particle we have reduced the 
minimum volume of fluid required for a characterisation of its vis- 
coelastic properties. Furthermore, our analysis relies on only a few 
key assumptions, particularly those relevant to the generalised 
Stokes-Einstein relation 1 , so should be robust within the realms of 
their validity. 

This technique will allow more detailed analysis of biological 
fluids, potentially for rapid and reliable diagnostics, and will increase 
the range of (bio)physical studies able to be carried out inside living 
cells 39 ' 40 and in chip-based devices. 

Methods 

Apparatus. The arrangement of optical devices is depicted in Fig. 1. The IR laser (IPG 
Photonics YLR-5-1064-LP, CW, 1064 nm with stable output at 2.0 W,P« 10 mWat 
the trap) is focused through an inverted high-NA microscope (Olympus, 1.3 NA, 
100 X objective and matched condenser lenses) to form the optical trap. A HeNe laser 
(JDS Uniphase 1137P, CW, 632.8 nm, P < 3 mW at focal plane) and plate are used 
for polarisation tracking. The AOMs (IntraAction DTD-274HA6, driven at 27 MHz 
with variable amplitude) allows the trapping power to be controlled. 

Fluids. Vaterite spheres were suspended in samples of Celluvisc (Allergan), a 10 g/L 
aqueous solution of CMC, after removal of the mother liquor through three repeated 
cycles of centrifugation, re- suspension in ethanol and sonication. 

Angle calibration. The transmitted HeNe light was passed through a polarising beam 
splitter cube and the resulting powers of both perpendicular beams measured. The 
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voltage difference generated (Ay) relates sinusoidally to the azimuthal angle of the 
vaterite probe by Ay oc sin 2(j). Introducing a quarter wave plate into the IR path 
causes the vaterite to spin. Recording the range of Ay allows an exact relationship 
between Ay and (f) to be found. 

Large scale rheometry. A parallel-plate strain controlled rheometer (ARES, 
Rheometric Scientific) with sample temperatures of approximately 8, 25 and 40°C 
achieved measurements of G* across the frequency range 1-200 rad/s by standard 
time-temperature superposition. 
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